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Abstract 

A linearized theory of the acoustics of porous elastic formations, such as rocks, saturated with 
two different viscous fluids is generalized to take into account a pressure discontinuity across the 
fluid boundaries. The latter can arise due to the surface tension of the membrane separating the 
fluids. We show that the frequency-dependent bulk modulus K(oj) for wave lengths longer than the 
characteristic structural dimensions of the fluid patches has a similar analytic behavior to the case 
of a vanishing membrane stiffness and depends on the same parameters of the fluid-distribution 
topology. The effect of the capillary stiffness can be accounted for by renormalizing the coefficients 
of the leading terms in the low-frequency limit of K (uj). 

PACS numbers: 43.20.Jr 
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I. INTRODUCTION 

Using the Biot theory [j| for the acoustics of fluid-saturated porous media as a starting 
point, it was possible |2| to formulate a simple theory for the frequency-dependent bulk 
modulus, K(u>), of a composite consisting of a homogeneous porous frame heterogeneously 
saturated with two different viscous Newtonian fluids. It is assumed that saturation occurs 
in "patches" of 100% saturation of one or the other of the fluids. In the limit when the two 
fluids are identical, the theory reduces to the low-frequency regime of the Biot theory. The 
mismatch of the local properties of the composite across the patch boundaries leads to a 
dispersion and attenuation mechanism jgj for acoustic waves, which is physically different 
from that of the Biot theory. In particular, the propagatory modes of the Biot theory, as 
applied to a fully monosaturated sample, are predicted to be nondispersive and nonatten- 
uative in the relevant regime, u <^ ub, where ujb is the Biot crossover frequency. In the 
case of patchy saturation, on the other hand, the mismatch of the elastic properties at the 
patch boundaries leads to a redistribution of the fluids between the patches when the sys- 
tem is perturbed. Because of the viscosity, this mechanism governs energy dissipation and, 
therefore, the attenuation and dispersion of the propagatory modes. There is a crossover 
frequency, u c , below which the pore pressure can equilibrate across the patches and above 
which it does not. 

This mechanism is particularly effective when there is a large contrast between the two 
fluids as occurs when one of them is a gas. In such a system, which was further studied in 
Ref. P| by analyzing measurements of Ref. Q, the gas saturated patches of the sample may 
often be approximated as vacuum. The relevant formulae become quite a bit simpler in this 
case. 

While the fluid redistribution across the patch boundaries was of a central importance 
in Ref. j^J, its theory failed to account for the membrane surface tension which can arise, 
for example, due to the pinning of the boundary edges inside the pores. The latter will 
result in the formation of a meniscus when the fluids are pushed between different patches, 
which in turn leads to a pressure drop between the patches due to the capillary forces. The 
present paper addresses this phenomenon. It is important to note that various additional 
capillary effects, such as the fractal nature of the fluid distribution in the pore space, can 
be important in realistic porous materials but lie outside the scope of this work. 
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Johnson |2j suggested that the patchy-saturation effect may be a dominant mechanism 
for the low-frequency (u <C wg) attenuation and dispersion in carbonate rocks, but can be 
dwarfed by microscopic "squirt" mechanisms in materials like sandstones. We indeed later 
showed U that the sonic and even ultrasonic measurements jj, 0, (| on various limestones 
can be successfully described by the patchy-saturation effect. In particular, we demonstrated 
there how one can extract information about the characteristic patch sizes and shapes in 
partially water saturated limestones by measuring the velocity and attenuation of acoustic 
waves. Extending our theory to take into account the capillary forces is necessary to improve 
our understanding of the acoustics in patchy-saturated materials. While the effect of the 
membrane tension at the boundary is but one of the important effects, investigating its 
effects can be a starting point for further studies of the capillary effects in porous media. 

The manuscript is structured as follows. In Sec. |H] we formulate the assumptions and 
limitations of the theory. Then, in Sec. IHII we in turn discuss the static, low-frequency, 
and high-frequency regimes of the theory in the general case of the patchy saturation by 
two arbitrary fluids. In Sec. IIVI an alternative formulation of the static and low-frequency 
limits is presented which is simple and physically transparent, albeit only valid for the case 
when one of the fluids is a gas. We compare our findings with exact numerical calculations 
in Sec. |V] for two simple geometries: that of the concentric spheres and the periodic slab 
geometries. Finally, in Sec. IVII a simple theoretical model for K(u) is presented enabling 
us to connect the low- and high-frequency asymptotes of the theory. Our results are then 
summarized in Sec. IVHI Two somewhat technical discussions, on the capillary correction to 
the high-frequency limit and analytic structure of K(uj), are given separately from the main 
text, in Appendices 1X1 and IB"! respectively. 



II. GENERALITIES 

As in Ref. Q], we consider a porous rock fully saturated with two different fluids, one 
having global saturation Si and the other S 2 = 1 — Si. Each point of the sample is saturated 
with one fluid or the other and the local elastic properties of the composite are assumed to 
be governed by the usual Biot parameters of the rock skeleton and the locally present fluid. 
The only exception is the patch and the sample boundaries which deserve special treatment 
(in the following). While we will talk about various limits of the frequency dependence of 
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the bulk modulus K(u), the Biot theory is always assumed to be in its low- frequency regime. 
The Biot crossover frequency 

= -j (1) 

(using the notation of Ref. |3|), defined in terms of the fluid density pj, the viscosity 77, the 
tortuosity of the pore space a^, the dc permeability k, and the porosity 4>, is thus the upper 
frequency limit of our theory, i.e., it is assumed that uj <C ujb- Biot equations of motion 
within each fluid patch |1| then simplify to 

V-r = 0, (2) 
Vp f = -£-(U-u), (3) 

where all quantities vary as u(r)e~ tuJt . Here, pf is the acoustic component of the pore 
pressure, U and u are the displacement of the fluid and solid phases, respectively, and r 
stands for the total stress due to both (solid and fluid) phases. Note that Eqs. (J2J), (J3J) do not 
contain the inertia terms, which is a relevant approximation as the attenuation/dispersion 
mechanism of our theory is governed by the slow diffusive mode. The pore pressure pf and 
the total stress r are related to the deformation strain in the solid phase, = (l/2)[uij + 
Uj ti \, and to that in the fluid phase, En = Uij, by jl| 

T = [(P + Q- 2N)e u + (R + Q)E U ] k 3 + 2Ne ij , (4) 
P f = —[Qeu + REu], (5) 

where the coefficients P, Q, and R depend on the bulk moduli of the solid phase, K s , the 
porous frame, Kb, and the pore fluid, Kf, as well as the porosity and the shear modulus of 
the porous skeleton, N (see, e.g, Ref. ^]). The fast compressional and shear velocities of the 
Biot theory are nondispersive and nonattenuative whereas the slow compressional wave is 
diffusive in the low-frequency limit considered here. Furthermore, we assume that the wave 
lengths of the fast and shear waves are larger than the characteristic fluid-patch dimension 
L, i.e., uj <C oj x , where 

2nV sh 

Wx = — — (6) 

with V^h being the (slower) shear-wave velocity. The latter assumption will allow us to 
define the effective bulk modulus K(uS) == —V(P e /5V), where (5V /V)e~ lu)t is the oscillatory 
fractional volume change in response to an external normal stress, r - n = —P e iie~ tuJt , n 



being the outward normal of the surface. Finally, the sample is considered to be isotropic 
and homogeneous, except for the fluid patches. 

Due to capillary forces, the fluid pore pressure is in general discontinuous across patch 
boundaries. Absent any applied stress, the pressure discontinuity across the boundary, i.e., 
the capillary pressure, — Pf\ is related to the mean curvature of the interface within 
each pore and to the surface tension via the Laplace- Young equation. In a completely 
different context than we are considering in this article, one is led to the equations of two- 
phase flow in which the partial saturation is generally treated as a well-defined function 
of the capillary pressure. See, e.g., Bear jjj. These concepts make no distinction between 
partial saturations that are patchy and those in which each pore may be partially saturated. 
Indeed, P^ and Pj 2 ^ are each generally considered to be macroscopic quantities defined at 
each material point throughout the system. 

In the present context, however, we assume the fluid patches are perturbed from equilib- 
rium by the application of an additional small- amplitude acoustic wave. The changes in the 
pore pressures are small compared to the equilibrium values: |p/(r)| <C \p9 {j = 1,2}, 
where p/(r) is the acoustic component of the pore pressure. On a microscopic level, the 
meniscus bends in and out, relative to its rest position. On a macroscopic level, there is an 
oscillatory displacement of fluid relative to the solid frame, at the boundary of the patch. 
This effect is conventionally captured by postulating the following boundary condition at 
the interface between two patches: 

pf - v f = W<j)(U - u) • n . (7) 

Here, pj is the acoustic pore pressure on the jth side of the boundary surface and n is an 
outward normal from patch 1 to patch 2. Un and un represent the components of the fluid 
and the solid displacements, respectively, that are normal to the patch surface; each of these 
entities is separately continuous across the patch boundary. Thus, Eq. (J7J) simply states that 
the additional amount by which the fluid interface bulges is proportional to the additional 
capillary pressure induced by the acoustic wave; the spring constant, W, is often called the 
"membrane stiffness". Indeed, Nagy and Blaho have shown that Eq. (JZ|) follows from 
the Laplace- Young equation applied to the case of pores that are circular cylinders. In that 
example they showed that W = cr/k, where a is the surface tension between the two fluids. 
This specific example illustrates the point that surface tension manifests itself as a stiffer 
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membrane in smaller pores than in larger ones. More generally, the membrane stiffness is 
predicted to take the form of 

W = sa/k , (8) 

where s is a shape factor which becomes much less than one as the pore structure becomes 
more irregular 11(1 . 

Our previa resets (aofiasoa Q, Tser— and Jofiusou Q) were specific to the ease 
W = 0. Here, we treat W as a real-valued, frequency-independent, phenomenological pa- 



rameter of the theory, to be measured independently in a given context. Nagy and Blaho 10] 
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have done just that for various porous media. See also Ref. jlll j. Furthermore, Nagy 
successfully used Eq. (J7J) to analyze his data on surface acoustic modes of fully-saturated 
porous media. In a somewhat different context, Eq. (JZj) is used by Liu and Johnson 
to describe the effect of an elastic mudcake on tube waves in permeable formations. The 
approximate validity of Eq. (j7J) is established on firm theoretical and experimental grounds. 
It is the essential purpose of the present article to investigate the effects of a nonzero value 
for W on the acoustic properties of patchy-saturated samples. 

Eq. (J7J) can be combined with Eq. (jSJ) to generalize the latter to hold at the interfaces as 
well as inside the patches: 

Vp/ + W r 0[n-(U-u)]fi(J(fl) = ^5(U-u), (9) 

where R(r) is the shortest distance from a given point r to the patch boundary and S(R) is 
the Dirac's 5 function. In addition, we assume the no net flow condition over the surface of 
the sample: 

JdS[V-u]-h = 0, (10) 

which is valid either for sealed outer boundary or for the periodic boundary conditions. (The 
integration is taken either over the surface of the sample, for the sealed-pore boundary con- 
dition, or over a unit cell, for the periodic boundary conditions.) Eq. (fTU|) can be converted 
to a volume integral using Gauss' theorem: 

(en) = (E u ) • (11) 

To summarize this section, the present theory supercedes that of Ref. |2j in that the 
assumption of vanishing capillary forces is now relaxed. The capillarity is accounted for by 
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considering the pressure drop across interfaces which is governed by membrane stiffness, see 
Eq. ( j7| l. The theory of Ref. |2| is thus recovered in the limit when W — ► 0. 

III. HIGH AND LOW FREQUENCY LIMITS 
A. Static regime 

In the static limit, it is possible to get an exact closed- form analytic expression for the 
compressive modulus, K , defined by Eqs. ©-Q- Hill jlj] has developed a well-known 
exact result for the bulk modulus of a (nonporous) composite in which the shear modulus is 
everywhere constant, though the bulk modulus varies from point to point. See also Ref. (15] ]. 
We extend that theory to the present case by looking for solutions of the form 

U{ OiXi + X,i i (12) 

Ui = axi + 4> ti . 

It is understood that \ as well as the normal displacement x,i n %i are continuous across the 
patch boundary; similarly for ip and ip^rii. Inasmuch as the first terms in Eqs. (fT2*|) are 
special cases of the second terms, one is free to choose any value for a, including zero. It is 
slightly convenient, however, to choose 

3a = (e„) = (En) ■ (13) 

Eq. (J3J) can now be written, within any patch, as 

Tijj = [{P + Q)x,n + (R + Q)tn], = • (14) 

This equation is satisfied if the Laplacians of x an d of ip are constant within each patch: 

. Oi re patch 1 
a 2 r G patch 2 

, 6i r G patch 1 

tn = T , • (16) 

02 r G patch 2 

The four constants, {ai, a2, &i, b 2 }, are determined by various ancillary conditions. First, the 
volumetric strains within the jth patch may be written as 

a) 

a (17) 
E}p = 3a + b,. 
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Therefore, Eq. (JT3*|) gives two conditions on the constants: 

Sjai + £202 = , 
Sjbi + S 2 b 2 = . 



(18) 



Next, there is the condition that the traction, r ■ n, is continuous across the boundary 
between patches. Hill jlj] has pointed out an identity satisfied by the discontinuity across 
the patch boundary of the second derivatives of \ which can be rewritten for our purposes 

as 

[xS-xSh = («i-"*)*> (19) 

where the superscript on x refers to the kth patch. Therefore, the discontinuity in the 
traction is 



T, 



(1) 



-(2) 



n ; 



{[{Pi + Qi~ 2N)(3a + oi) + (P x + Q x )(3a + &i)] 
- [(Pa + Q2- 2N){3a + o 2 ) + (P 2 + Q 2 )(3a + b 2 )} 
+ 2N(a 1 -a 2 )}n i . 



(20) 



(Pk,Qk,Rk refer to P,Q,R evaluated with respect to the saturating fluid in patch k.) In 
order for Eq. (J2*U|) to vanish across the patch interface, it is necessary that the quantity in 
curly brackets vanishes: 

(Pi + Qi)(3a + 01) + (Pi + Qi)(3a + 61) = (P 2 + Q 2 )(3a + a 2 ) + (P 2 + Q 2 )(3a + b 2 ) . (21) 

This, then, is the third condition on {01,02,61,62}. 

The fourth and final condition follows from Eq. (J7j). The pore pressure within the fcth 
patch is given by Eqs. (0) and (fTTj) : 



(fc) 
Pf 



i[P fc (3a + b k ) + Qfc(3a + a k )} . 



(22) 



We see that the pore pressure is constant within each patch, as required by Eq. Q in the 
static limit. The left-hand side of the boundary condition is constant over the interface 
between the two patches. Consider the case of a finite volume with sealed boundaries. 
Integrating Eq. (0) over the surface between the two patches, one finds 



' (1) (2) 
Pf Pf 



A = W<p dA(U-u) -n, 

J A 



(23) 
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where A is the area of the interface between the patches. This integration can be extended 
over the entire bounding surface of fluid 1, say, because there is no flow across the outer 
boundary: 

' ' A = W(j) [ dA(U - u) • n . (24) 



' (1) (2) 
Pf Pf 



Now one may use Gauss' theorem: 



Pf 



(2) 
Pf 



A = W(j) / dV(E u - e u ) = W(j)S x V 

JVi 



(25) 



The volumetric strains and the pore pressure can be related to the unknown constants by 
means of Eqs. (JT7J) and (|2*2*|). respectively. The result for Eq. (J23j) is then 

fli(3a + h) + Qi(3a + a x ) - [R 2 {3a + b 2 ) + Q 2 (3a + a 2 )] = Wfa - b 1 )S 1 , (26) 

where 

W = W(p 2 j (27) 

in terms of the sample volume, V, and the bounding area between the patches, A. If the 
derivation is repeated for the case of periodic boundary conditions, Eq. (j2EJ) still holds; in 
that case, V is the volume of a unit cell and A is the interface area between the patches 
within the unit cell. In either case, if the derivation is repeated by closing the surface over 
fluid 2, instead of fluid 1, the right-hand side of Eq. (|2*S|) is replaced by —W(a2 — &2)5 , 2, 
which is equivalent via Eqs. (JTHj). 

The four linear equations, (fTHj) . (|2~T|) . and (J2HJ), can be inverted analytically to get the 
ratios {ax, a 2 , &i, &2}/a;; in practice, this is more easily done numerically. The effective static 
modulus is simply related to these constants. Following Ref. we identify the external 
pressure, P e , with the average compressive stress in the system: 



(2) 



(28) 



where the compressive stress in the kth patch follows from Eq. (fl|): 



t, 



(k) 



^P k + Qk--N^ja k + (R k + Q k )b k + 3K^a . (29) 

Here, K^q = P k + 2Q^ + — (4/3) AT is the Biot-Gassmann modulus for the kth. patch. 
The effective static modulus for the system, K = —P e / (en), is thus given by 



K 



I* 



+ S 2 



Px + Qx- \n] ^ + (R t + go - + 3^ 

6 J a a 



P2 + Q2-^n)^ + (r 2 + q 2 )^ 

6 J a a 



3K 



(2) 
BG 



(30) 
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Within the context of the model, Eq. (|5U|) is exact. It depends upon the usual Biot 
parameters, the saturation values, the membrane stiffness, W, and it also depends upon the 
surface to volume ratio of the patch interface, A/V. There are three interesting limiting 
cases of Eq. (130)) : (1) If the two fluids have identical bulk moduli, then 

J m K „ K ° = K ™ = K ™- < 31 > 

K f ~ K i 

independent of the value of the membrane stiffness parameter, W. When the sample is 
compressed, the pore pressure rises equally in the two patches with the result that there 
is no relative flow between fluid and solid. (2) If the membrane stiffness is large enough, 
then each fluid is locked within its own patch and the static modulus is equal to the Biot- 
Gassmann-Hill result: 

lim K = ^bgh , (32) 

where 

1 = Sl + S2 (33) 

^bgh + (4/3)7V K^ + (4/3)N + (4/3)N ' 1 } 

(3) If the membrane stiffness is zero, then the static modulus is given by the Biot-Gassmann- 

Woods result as was proved previously j^: 

Urn K = K BGW , (34) 

where -Kbgw = K~BG(Kf S ) and Kj s is an effective fluid modulus given by the Wood's law 
expression: 

' *' +t!t ■ (35) 



Kf Kf Kf 

In Fig. show the behavior of K as a function of W. The parameter set [16] is based 
on that of Ref. [2j. The sample is saturated with water, S w = 0.7, and with air, S a = 0.3. 
If the membrane stiffness is small, then Eq. ()34j) holds whilst if it is large, Eq. ()32j) is valid. 
The crossover occurs when W is approximately equal to the modulus of the stiffer fluid, 
water. In Fig. |2Jwe show the saturation dependence of the static modulus, K , as well as 
that of i^Bcw an d i^BGH- Even when the membrane-stiffness parameter, W, is independent 
of the saturation, the parameter W will in general have a saturation dependence reflecting 
the change in the interfacial area. Nonetheless, for this plot, W is held constant at a value 
W = 1 GPa. 
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B. Low-frequency limit 



Here we wish to find the behavior of K(u) as it approaches Kq. Using the method devel- 
oped in Ref. J^j, we consider two equivalent expressions for the average energy dissipation 
per cycle. The macroscopic expression for the dissipation power averaged over one cycle is 



1 r du* 
-Real / dS—— 

2 J dt 



t ■ n 



-Real 
2 



\Pf 



(36) 



where the integral is taken over the bounding surface of the sample. The microscopic 
expression is 

1 f (d\J* du*\ 
P = --Real / dV0 I— - — ) ■ Vp, . (37) 



dt 



In Eq. ()37j) the integration is taken over the volumes of each patch, specifically excluding the 
bounding surfaces of the patches; Vp/ has a 5-function contribution on the patch surface, see 
Eq. One may then substitute Eq. (j^J into Eq. ()37)1 to arrive at the following microscopic 
expression for the dissipation: 



P 



--Real / dV<jriu{\5* - u* 



iujipr] 
k 



(U 



u 



2k 



J dVr](r) |U(r) - u(r) 



(38) 



In this expression, the integrand does not contain any 5-function contributions. It is imme- 
diately clear that, as expected, the low-frequency behavior of P is quadratic in frequency. 
Therefore, the low-frequency expansion of K{u) is given by 



lim K(u) = K 1 - icoT + 0(icu) 



(39) 



where the parameter T is determined by substituting Eq. into Eq. and equating it 
to Eq. flUED: 



T 



dVr](r) 



U(o) ( 



u (o) (r] 



(40) 



kV\P e \ 2 

which is expressed in terms of the static displacements, U^ -* and u^. 

We know from Eqs. (|12|) that the displacements in the static limit are given by the 

gradient of certain functions. Therefore, we may define an auxiliary function, $, by 

Tj(o) _ u (0) 



V$ = — rj(r) 



(41) 



Within each patch $ = i][x — ip]/ P e + const. The function $ is continuous across a patch 
boundary (otherwise there would be a 5-function contribution to the displacements); the 
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normal component of the vector — V$/?](r) is also continuous; and $ is either periodic or 
V$ = (in the case of the sealed-pore boundary condition) on the outer boundary. The 
function $ satisfies the differential equation 

This equation is analogous to an electrical conductivity problem in which <&(r) plays the 
role of a potential, rj (r) is the local resistivity, and g(r) is the distributed current source. 
Electrical neutrality is guaranteed by Eq. (jllj) : / dVg(r) = 0. 
In terms of $, Eq. (|4U|) may be written as 

r ^/ d %M< v *> 2 < (43) 

which may be integrated by parts: 

(The surface term vanishes identically for either the sealed or the periodic boundary condi- 
tion.) The source function, g(r), is constant within each patch but in general is discontinuous 
across the patch interfaces. From Eqs. (|T2|) . and (JT7|) one has 

(6i-0!)/P e rG patch 1 
0(r) = <{ . (45) 

(b 2 — a 2 )/P e r ^ patch 2 

Thus, the determination of the parameter T which governs the low-frequency approach to 
the static limit, Eq. ()39|). requires solving the differential equation in which the piecewise 
constant sources, Eq. (j45jl . are known by virtue of the solution for the static regime, and 
imposing one of our two no net-flow boundary conditions. This situation is only slightl 
different from that encountered when the surface-membrane effect is ignored altogether 
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C. High-frequency limit 

Here we consider frequencies that are high in the sense that the pore pressure does not 
have time to equilibrate from patch to patch but we are still in the regime where to ujb 
and u uj x . In this high frequency limit, to —>• oo, the fluids within each patch are locked 
with the solid frame, e« = En. This conclusion is independent of the interface-stiffness 
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parameter, W. We have a situation which satisfies the assumptions of Ref. jl^ : The shear 
modulus is a global constant and the bulk modulus is constant within each patch. Hence, 

Koo d ^ f lim K(u) = K BG}1 (46) 

for any value of W, where K BG n is given by Eq. Recall that we found the same result 

for the static compressive modulus, Kq, when the patch boundary is locked by the large 
capillary tension, W — > oo, see Eq. 

For a large but finite frequency, the correction to is provided by the diffusive modes 
at the patch boundary [2] . In terms of a local coordinate x normal to the interface, the pore 
pressure near the boundary, x = 0, is given by 

| p fl + A x e- iq ^ x , x < (patch 1) 
[ Pf2 + A 2 e iq2X , x>0 (patch 2) 
Here, pfj. is the pore pressure inside the kth patch in the limit of infinite frequency and 



q k = \JiujD k (48) 

is the slow-mode wave vector in terms of the diffusion coefficient see Appendix EI I n 
the case of a finite membrane stiffness, Pf(r) is in general discontinuous across the patch 
boundary. We can, nevertheless, show (see Appendix EJ that the pore pressure discontinuity 
8pf = f \p^p —pf\ = 0(uj~ 1 I 2 ) Apf = f \pfi — pf 2 \ = 0(uj°), for large frequencies, and thus 
the leading correction to can be found by setting 5pf = 0, i.e., assuming continuity of 
the pressure, Eq. P7j). Requiring, in addition, continuity of the fluid flux, we can then find 
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coefficients A 12 in Eq. (jjTj) and reproduce the result of Ref. 

lim K(u) = K BGH \l - {iuo)- l/2 G + ■■■ , (49) 
where 

G — fe^BGH /Wj (5Q) 

in terms of the would-be drop in the pore pressure, Apf, relative to the applied external 
stress, P e : 

Ap f = (R 2 + Q 2 )[K 1 + (4/3)iV] - (gi + Q 1 )[K 2 + (4/3)7V] 
P e <f)S 1 K 1 [K 2 + (4/3)N}+(j)S 2 K 2 [K 1 + {4/3)N] ' 1 j 

Notice that the membrane-stiffness parameter does not enter any of Eqs. (|49|) - (j51|) . so that 

the capillary forces not only have vanishing contribution to the high-frequency compressive 

modulus (j4T?j) but also to the leading correction to it, Eq. (TE)j) . 
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IV. LIQUID AND VACUUM 



In Sec. IHII we have seen that whereas the two terms in the high frequency limit, Eqs. 
are identical to those derived earlier when W — 0, the two terms which describe the low 
frequency limit, Eqs. (|39|). are now considerably more complicated. The prescription for 
finding the static bulk modulus, K [Eq. (JHUj)]. and the coefficient of the leading low-frequency 
correction, T [Eq. (|44|) ]. follows from the solution of the system of equations for parameters 
{ai, 02, b±, 62}; no simple analytic form is deduced and no transparent physical interpretation 
is readily available. We therefore offer an alternative derivation, which is both simple and 
physically appealing, albeit valid only for the case when one of the two fluids is a gas taken 
here to be the vacuum, and the other is a liquid whose properties are similar to water or oil. 
(A gas at standard temperature and pressure is four orders of magnitude more compressible 
than water or oil and so the approximation is a good one except for gas saturations less than 
~ 10- 4 .) 

As in Sec. IIII| the starting point is Hill's theory If we denote the bulk modulus of 
the patch saturated with the liquid by K', then the static modulus of the composite is given 
by 

1 S 1-S 

+ 1? TT77ov?\7 . (52) 



A^ + (4/3)A K'+(4/3)N K b +(A/3)N 
where S is the fluid saturation and Kb is the bulk modulus of the dry frame, as before. 
K' can be found by extending the Biot theory: Since the liquid does not bear any shear 
stresses, the static compressive modulus of the liquid-saturated patch must be given by the 
Biot-Gassmann formula, 

K' = K BG (K' f ) , (53) 

for some effective fluid bulk modulus, K^. This effective modulus characterizes the "effort 
of squeezing" the liquid out of the pore space. In particular, two limiting cases are obvious. 
(1) K'r vanishes with the membrane stiffness: 

lim K'r = (54) 
w->o J v ' 

and (2) it reduces to the bare fluid modulus for infinitely large stiffness parameter: 

lim K' f = K f . (55) 
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If by inducing the pore pressure pf, the fluid-filled pore-space volume is reduced by SV, then 
K'j is defined by 

If the fluid volume squeezed out of the pore space into menisci formed at the patch surface 
is SV', then 

Pf A = WSV , (57) 
which follows from Eq. (|24|). Finally, we have 

for the bare bulk modulus of the pore fluid. Combining Eqs. ()56|) - fl58|) we then arrive at 

11 A 
K r f ~K~ f + W0SV ' (59) 

which satisfies Eqs. ()54j) and (|55|) in the two limits given by Eqs. (|54|) and ()55jl . As a 
consequence, in the limit W — > 0, K' = K b (which is a specific case of a general result holding 
for arbitrary two fluids: K' = K BGW ), and for W — > oo, K' = Kbgh- Equations (jB^j) . (fo^j). 
and (J59)) complete the derivation for the static limit. This result for the static modulus is 
also plotted in Figs. Q and El where it is seen to be identical to that computed from the full 
formalism. 

It is as straightforward to find the parameter T governing the low-frequency behavior. T 
is given by Eq. (}4"4"|) after solving Eq. (|42jl . when the function g(r) is known. In the fluid-filled 
patch, the volumetric strain of the solid can be found using Hill's theory once again (see 
also Eq. (35) of Ref. U): 

[K b + (4/3)iV]P e 

11 K'K b + (A/3)N[SK' + (l-S)K b ]' 1 ' 

Using Eq. © and 

Pf A = W(j)SV{E u - e u ) , (61) 
provided by Eq. (|25|) . we then obtain 

Q + R 

Eu — eu = — = €u , 62 

WS + R 

which together with Eq. (|6T)J) defines the value of g, Eq. (jlSJ), inside the fluid patch: 

a = ( Q + R \ [Kb + (4/3)iV] f63^ 

y \WS + RJ K>K b +(4/3)N[SK>+(l-S)K b )- 1 ; 
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Here, quantities Q and R are evaluated using the bare fluid bulk modulus, Kf. Defining an 
auxiliary function $ in the space occupied by the fluid patch, Vf, by 

V 2 $ = -1 r 6 V f (64) 

subject to the patch boundary condition $ = as well as the usual no net-flow condition 
on the external surface of the sample, we find for the parameter T, Eq. (|44j) . 

m 

where the coefficient has the dimensions of a diffusion constant and is given by 

DT = lWiuf ' (66) 

The Poisson size of the fluid patch, If, is defined by 

S-VjL,**- (67) 

The sole effect of the membrane stiffness, W, in this limit, is thus to renormalize the coeffi- 
cient Dt from that which held when W = 0. We have verified numerically that the results 
derived in this section are identical to those presented in Sec. IIHI when one of the fluids is 
taken to be a vacuum. 



V. EXAMPLES 

It is straightforward to repeat the calculations of Ref. |2| for K(uS) in which the patchy 
geometries are either that of periodic slabs or concentric spheres. In the slab geometry, region 
1 is a layer of thickness L\ and region 2 is a layer of thickness L2, periodically repeated. 
Here, Si = L 1 /(L 1 + L 2 ) and A/V = 2/(L 1 +L 2 )- m the concentric-spheres geometry, region 
1 is a sphere of radius R a surrounded by region 2 of outer radius R},: Si = {R a f Rb) 3 and 

A/v = ml/Rl 

Equations (B.4) and (B.8) in Ref. |2j are still correct as written for the slab and the 
sphere geometries, respectively. However, the pore pressure is no longer continuous across 
the patch boundary but it obeys Eq. instead, leading to a slightly different matrix 
equation to be solved. The results of such numerical calculations are presented in Fig. 01 
using the parameter set of [r| at a water saturation S w = 0.9. The horizontal dashed lines in 
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the left column are Kbgb., K q , and Kbgw from high to low, respectively. At low frequencies, 
we see that the numerically exact results do indeed approach K Q , as per Eq. Abgw 
has the same value in all three geometries, because the saturation values were chosen to be 
the same. Similarly for A"bgh- The values for K , however, are different in the three cases 
because the surface to volume ratios for the different patch geometries are different. The 
parameter W in Eq. ()26|) has a different value in each case and so, therefore, do the values 
of K . 

As we mentioned earlier, the low-frequency parameter T is determined by a differential 
equation (|42jl . which is, in form, identical to that which existed when W = Q]. This 
means that the previously derived results for the slab and the sphere geometries are virtually 
unchanged. For the sphere, we have from Eq. (40) of Ref. 



3 d2 i r _ ro., v -, „ ( o v , i „ \ „ 1 d2 ]-)'■'> o „ Jl j j 5 1 



-15 V2 g 2 (g 2 - gx)R 6 a Rl + 5g 2 [3rj 2 g 2 - (2 V2 + m)gi]RiRt - ^ 2 g 2 2 Rl} . (68) 



{Vi9i L l + Zr)ig x g 2 LlL 2 + ^ 2 g x g 2 L x L\ + J] 2 g\L^ . (69) 



Similarly, the expression appropriate to the slab geometry is Eq. (41) of Ref 

K n 6 2 



Qk{L x + L 2 ) 

Of course, the correct value for the static modulus, K , Eq. (JHUj) . and the correct values 
for gj, Eq. (|45|). must be used. Subject to this clarification, the low-frequency attenuation 
implied by the second term in Eq. (ffilj) is shown as a dotted line; we see that the numerically 
computed results asymptote correctly in this limit. It was noted in Sec. IHII that the high- 
frequency asymptote for K(u) is unaffected by a nonzero value for W. This limit is also 
plotted as dotted lines in Fig. El 

VI. THEORETICAL MODEL FOR K(u) 

A key observation by Johnson [2] was that all singularities and branch cuts of the bulk 
modulus K(oS) must lie on the negative imaginary axis, when analytic continuation into 
the complex plane is performed. This stems from the diffusive nature of the relaxation 
mechanism. In Appendix [BJ we show that this conclusion stays valid for a nonzero value for 
the membrane stiffness, leading us to suggest that the simple analytical formula proposed 
earlier [2] may also work in the present context. We wish to have a function which has 
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singularities only on the negative imaginary axis and which has the limiting properties of 
Eqs. (J3HJ) and (|39~|l . That is, we propose the following function should work well as compared 
against exact numerical results: 

K{u) = K OQ K oo~K 

i - c (i - V 1 - wc 2 ) 

where consistency with Eqs. and flT^ requires 



(71) 



and 



= (Koo - K,f = (K^ - go) r 
^ 2KqK 2 00 TG 2 2K T' { ) 

In Eq. f!70|) the branch cut in the definition of the square root function, \fZ, is taken to be 

along the negative real z axis. 

For the specific cases of the slab and concentric spheres geometries, the parameter T is 

computed directly from Eqs. (JUH|) . (joT?|) as the case may be. G is evaluated using Eq. (JoTty . 

The resulting predictions from Eq. (J7U|) are shown as dashed curves in Fig. (jSJ), where one 

can see that it does an excellent job of interpolating throughout the entire frequency range. 



VII. CONCLUSIONS 

In summary, we have extended the theory of Johnson to consider effects of the capillary 
forces on the acoustics of porous media saturated with two different fluids. Using a well- 
known expression for the pore-pressure discontinuity across the patch boundary, Eq. ((7|), we 
have shown that the finite membrane stiffness, W, does not affect the analytic structure of 
the frequency-dependent bulk modulus K{u). In addition, the high-frequency asymptote of 
K{ui) does not depend on W, as long as W is sufficiently small [see Eqs. flAlOj) . (jAllj) ]. so 
that the high-frequency limit of our theory can be reached at all. The effect of the capillary 
forces is thus to rescale the static bulk modulus, Kq, and the low-frequency correction 
parameter, T, Eqs. (|3U|) and (|44|) . respectively. Furthermore, there are no new geometrical 
parameters introduced by the membrane stiffness: For example, when one of the fluids is a 
gas, the two relevant topological parameters are the sample-volume to patch-surface ratio 
and the Poisson size of the fluid patch, //, Eq. (JoTj) . In the limit W — > 0, results of Ref. P] 
are recovered. In the limit W — > oo, K(u) = /Tech, Eq. (j3HJ), for all frequencies. 
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APPENDIX A: CAPILLARY CORRECTION IN THE HIGH-FREQUENCY 
LIMIT 

In this Appendix we show that the pore pressure discontinuity Spf across the patch 
boundary, see Eq. ([47)1 . scales as uj~ 1 ^ 2 at high frequencies, u> — > oo, so that the membrane 
stiffness does not modify the high-frequency limit. The pressure equilibration, Eq. (J47|) . is 
provided by the slow longitudinal waves on both sides of the patch boundary. In order to 
find equations of motion for the slow wave inside either patch, we take divergence of Eqs. (J2J), 
(j3J) and use relations ©: 



{P k + Qk)4 ] + {R k + Qk)E-, 



(k) 
u 



V 2 



Looking for solutions of the form 



Qk4l ) + RkE n 



(k) 



0, 

iu(f) 2 r) k 



.(*) 



E. 



(k) 



(Al) 
(A2) 



E 



(k) 



c ( fc ) e -i(q k x+uit) 



(A3) 



we obtain from Eq. (jAl|l the ratio of the amplitudes of the solid and fluid phases 

_ Rk + Qk 

Pk + Qk 

and then from Eq. (JA2j) the diffusion coefficient 

def iu k 



(k) 

c 2 



(A4) 



PkRk — Qk 



4 



(A5) 



f Vk P k + 2Q k + R k ' 

The above equations are valid on both sides of the patch boundary. The coefficients of 



pressure correction Ak, Eqs. (J47|) . at the interface are then given by 



A, 



Qkc[ k ^ + Rkci 



(fc) 



(A6) 



where we used Eq. (JSJ) for the additional pressure exerted by the slow waves. The pressure 
drop 5pf across the interface is related to the relative fluid displacement through Eq. ((7|), 
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and using Eq. (jSJ) we thus obtain 

5p f = (p fl + A 1 )-{p f2 + A 2 ) = -^-^A 1 . (A7) 

7]iUJ 

Finally, we require the fluid-flux continuity across the patch boundary: 

^4 1 = -*U 2 . (A8) 

Eqs. (jA4|) . ()A6|) - ()A8|) now constitute a linear system of 6 equations with six unknowns: c[ k \ 
c 2 , and Ak- The resulting expression for 5pj, Eq. (|A7j) . simplifies considerably if one patch 
is filled with a vacuum. In this case, the pressure discontinuity at the fluid boundary, relative 
to the difference in pore pressure relatively far from the interface is given by 

*EL = I = 1 (A9) 

Ap f 1 + rju/Wkq i + ^/i^ ' 

where the crossover frequency 



<*w = ^ (A10) 

is defined in terms of quantities corresponding to the fluid patch. When uj 3> ujw, 8pf oc 
uj^ 1 ! 2 and the pressure drop across the interface, Spf, is negligible in comparison with the 
pressure difference inside the patches, Apj, Eq. (JoTj) . 

For consistency of our theory, we require that this limit can be reached while we are 
still in the low-frequency regime of the Biot theory and the propagatory modes have wave 
lengths longer than the typical patch dimensions, i.e., we need 

lo w < (ub,u x ) . (All) 

Eqs. (|A10|) and (|A11|) imply that our high-frequency limit is valid only for sufficiently small 
membrane-stiffness parameter W. With the parameter set used to generate figure El we have 
uw/CZn) = 4.5 Hz. However, our theory makes no assumption about the value of u>w as 
compared against the diffusive crossover frequency and, indeed, the theory works quite well 
even if W is so large that u>\y is greater than the crossover implied by the intersection of the 
high and low frequency asymptotes. 

APPENDIX B: ANALYTIC STRUCTURE OF K(u) 

Here we generalize Appendix C of Ref. j^] to take into account the additional term due 
to the finite membrane stiffness W in Eq. Q. We want to investigate analytic properties of 

21 



W 2 k 2 



the real-valued causal response function K(t) in the frequency domain: 

K(u) = [°° K(t)e luJt dt. (Bl) 
Jo 

It is clear that K(u) is analytic everywhere in the upper-half complex uj plane, and we will 
further show that all the zeros, singularities, and branch cuts lie on the negative imaginary 
axis. 

In the following we consider an eigenvalue problem in which all quantities vary as 
u(r)e~ luJnt . Assume 

J dSu-T* • n = , (B2) 

as would be in the case of a zero, singularity, or branch cut of K(u) 2]. We can then show 
using Eq. (J2J) and the Biot relations between stresses and strains, Eqs. (@J) and ©, that 
/ dVpfTu is real- valued. This is true since all these equations are unaffected by the capillary 
forces, and we thus reproduce Eq. (C5) of Ref. 

1 



JdVp f 



rT, = dV 



-M 2 + ^(DD*)„ 

3 2N y 111 . 



(B3) 



K b /K s - 1 

where = — (l/3)rjj5y is the deviatoric part of the stress tensor. Finally, multiplying 
Eq. (JOJ) by U* — u* and integrating it over the entire volume of the sample, we obtain 



iuj. 



J d y^|u -u\ 2 = J dVp f {e* u - E*) + W<j> J dV\h • (U - u)\ 2 6(R) . (B4) 



As the volumetric strains, en and En, are given by linear combinations of pf and Tu [via 
inverting Eqs. (jlj), (jSJ], the right-hand side of Eq. (|B4|) is real- valued and, therefore, u n is 
pure imaginary. This completes the proof. 
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FIG. 1: Dependence of the static modulus on W for a sample saturated with water, S w = 0.7 and 
air, S a = 0.3. Parameters are listed in Q]. The two results for the static modulus, Eqs. (|3T)|) and 
(|52|). are seen to be nearly identical in this case in which one of the fluids is gas. -Kbgh > ^bgw- 

FIG. 2: The static moduli as a function of the water saturation, S w . Parameter set as in Fig. ^ 
W = 1 GPa is held constant. Again, Kq is calculated in two different ways. -Kbgh > -^bgw- 

FIG. 3: Dispersion (Realfi^u;)]) and attenuation [1/Qk = —Im&gK (u>) /~R.e&lK (uj)] due to the 
patchy-saturation effect for an air/water combination, S a = 10%, in three different situations. 
Solid curves are numerically exact solutions, dotted curves are the high- and low-frequency limiting 
expressions, and dashed curves represent the analytic expression, Eq. (|70|) . Horizontal dashed lines 
are (in order of decreasing value) Kbgh, Kq(W = 30 GPa/m), and Kq(W = 0). (A) Top Row: The 
gas is located in an inner sphere, R a = 4.642 cm, surrounded by a shell of water, = 10 cm. (B) 
Middle Row: The roles are reversed, R a = 9.655 cm, R^ = 10 cm. (C) Bottom Row: A p eriodic 
slab geometry, L a = 2 cm, L w = 18 cm. The other material parameters are listed in [16j. In all 
cases, the membrane stiffness is W = 30 GPa/m. The values of Eq. (|72j) . are shown. 
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Membrane Stiffness Effect 
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Figure 1. Y. Tserkovnyak and D. L. Johnson, J. Acoust. Soc. Am. 
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Saturation Dependence of Moduli 
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Figure 2. Y. Tserkovnyak and D. L. Johnson, J. Acoust. Soc. Am. 
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Figure 3. Y. Tserkovnyak and D. L. Johnson, J. Acoust. Soc. Am. 
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